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ABSTRACT 


The  small  perturbation,  two-dimensional  transonic  equation 
is  manipulated  with  a  separation- of -variables  approach  to 
obtain  two  ordinary,  nonlinear,  differential  equations. 
Numerical  integration  of  these  differential  equations  results 
in  new  transonic  boundary  surfaces  for  planar  external  flows. 
A  key  ingredient  in  these  solutions  is  the  identification  of 
dependence  of  two  integrations  constants,  a  and  /?,  on  the 
parameter  (1-M^^)  .  The  anticipated  behavior  for  both  the  Mach 
number  and  the  pressure  coefficient  is  used  as  a  guide  in  the 
actual  selection  of  the  adjustable  constants  in  the  problem. 
The  physical  reality  of  our  boundary  surfaces  is  examined  by 
displaying  the  boundary  conditions  they  satisfy.  The  strictly 
sonic  flow  (Mo„=1.0)  has  an  analytic  representation 
corresponding  to  a  divergent  surface  which  goes  supersonic. 
This  sonic  solution  is  compared  with  an  Euler-CFD  approach 
confirming  the  validity  of  our  results  over  the  region  where 
small  perturbations  apply.  Solutions  are  also  shown  for 
Mo„=0.8, 0.9, 1.1,  and  1.2  .  These  results  are  consistent  with 
known  behavior  for  both  subsonic  and  supersonic  external  flow. 
Since  the  results  of  this  work  yield  actual  transonic 
contours,  we  can  examine  shockless  surfaces  for  design 
applications.  The  possibility  of  starting  with  transonic 
surface  is  of  interest  to  present  day  CFD  approach.  Finally  an 
entire  transonic  upper  surface  is  presented  for  Meo=0.8,  by 
patching  a  subsonic  Mach  number,  which  reaches  a  plateau  at 
M=1.0,  with  a  sonic  flow.  This  patching  requires  the  careful 
interpretatio7i  of  a  nondimens ional  reference  length,  called 
Yo,  which  is  a  function  of  M„. 
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I.  INTRODUCTION 

The  governing  nonlinear  equations,  together  with  realistic 
boundary  conditions  have  always  been  the  main  difficulty  in 
transonic  aerodynamic  analysis,  [Ref.  1]  .  While  the  hodograph 
transformations  approach  linearizes  these  equations,  the 
transformations  are  confined  to  a  limited  number  of  transonic 
flows  due  to  complicated  boundary  conditions.  The  small 
perturbation  approach  simplifies  the  full  potential 
Equations;  furthermore,  like  other  predictive  methods,  it 
succeeds  in  retaining  the  great  sensitivity  of  transonic  flows 
within  small  perturbations.  This  paper  expands  on  the  exact 
solution  of  the  small  perturbation,  nonlinear,  two 
dimensional,  transonic  equation,  using  the  guidelines  and 
tools  given  by  Biblarz  0.,  [Refs.  2  and  3].  Starting  v/ith  a 
separation- of -variables  approach.  Chap.  Ill  reveals  two 
nonlinear,  ordinary,  differential  equations  which  lead  into 
tvjo  exact  solutions.  Chaps.  IV  and  V  describe  the  numerical 
integration  of  these  implicit  solutions,  which  finally  yield 
the  boundary  surfaces.  Later,  these  boundary  surfaces  are 
found  to  satisfy  the  boundary  conditions  for  a  two  dimensional 
surfaces.  These  findings  are  then  compared  with  a  finite 
difference  solution  of  the  Euler  Equations  for  the  sonic 
solution.  The  usefulness  of  the  shockless,  transonic,  boundary 
surfaces  appears  once  transformed  to  dimensional  body 
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surfaces;  "patching"  (i.e.,  translating),  the  dimensional, 
subsonic,  body  surface  of  Cp=0  at  the  upstream  inflow,  with 
the  dimensional,  sonic  body  surface,  enables  us  to  form 
complete  transonic  upper  surfaces  of  interest  for  design. 
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II.  THE  TRANSONIC  EQUATION 


The  small  perturbation  nonlinear,  two  dimensional 
Transonic  Equation 

+$  =  $  (1) 

V  J.  'ifyy  ^ 


is  described  by  the  velocity  potentials 

U  =  V  =  <^y  (2) 


Multiplying  eqn.  1  by  1/U„  enables  us  to  transform  the 
Transonic  Equation  to 

where  the  nondimensional  velocity  potentials  are 


4)x  = 


_u_ 


(4) 


The  Meyer  solution  for  the  de-Laval  nozzle  defines  a 
velocity  potential 


(})(^/y) 


=  b 


/Jt  tU  u  4 


(5) 


which  satisfies  the  nondimensional  form  of  equation  (3)  for 
M„=1.0  .  The  Meyer  solution  belongs  to  a  general  class  of 
sonic  flows  as  reported  by  Guderley  and  Yoshihara,  [Ref.  4], 
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The  transonic  equation  (1)  may  be  modified  to  a  more 
useful  form  if  multiplied  by  Mj^{y+1)  /U„,  instead  of  just  1/U- 


,,  to  give 


+  <^yy  = 


(6) 


Here  the  modified  velocity  potentials  are 


4)^  =  MUy^l)^ 


u 


i^y  =  M^(Y+1)-^ 


(7) 


For  a  given  problem,  and  y  are  constants,  and  eqn.  6  is 
really  less  cluttered.  However,  notice  that  M„^(7+l)/U„  is 
already  given  in  the  derivatives  of  4>  in  eqns.  7,  and 
therefore  appears  in  the  expressions  of  Cp  and  the  boundary 
surfaces  as  will  be  shown  later.  This  form  of  the  modified 
transonic  equation  (6)  will  be  examined  in  the  next  chapters. 
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III.  THE  GOVERNING  DIFFERENTIAL  EQUATIONS 


An  exact  solution  to  the  modified  transonic  equation  6  has 
been  given  by  Biblarz,  [Ref.  2] ,  starting  from 

4)(x,y)  =  (t)^(x,y)  +  {1-mI)  X  (8) 

where  0®(x,y)  refers  to  the  sonic,  small  perturbation  solution 
(M„=1.0) .  In  other  words,  "any  flow  which  satisfies  the  sonic 
equation  also  satisfies  the  full  transonic  range  with  similar 
boundary  conditions,"  [Ref.  2]. 

An  exact  solution  may  be  obtained  by  using  the  separation- 
of -variables  approach  with  the  potential  function  4){x,y)  of 
the  form 

4)(^/y)  =  ^ix)i\{y)  +  [I-mI)  X  (9) 

Substituting  the  above  ^  function  in  the  modified 
transonic  equation  6  results  in  two  ordinary,  second  order, 
nonlinear  differential  equations 

^  =  0  (10) 
dx  dx^ 

=  0  (11) 

dy2 

where  X  is  a  separation  constant. 
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A  simple  solution  to  the  first  differential  eqn.  10  is 
obtained  by  multiplying  both  sides  by  d|/dx. 


a?)  =0 

dxydx  dx^  dx 


or 


_d 

dx 


3\dxj  2^ 


=  0 


Thus 


dx  \  2 


dx  = 


d^ 


M  —V  +  a 

N  2 


and 


X 


-  ^0  =  /- 


d? 


3A 


$2  +  cj 


where  a  and  are  integration  constants. 

A  similar  procedure  can  be  performed  to  solve  the 
differential  eqn.  11  yielding 


dy\  dy^ ) 


- 

dy 


=  0 


or 


d  '  1/  dr] 
dy[  2  \  dy) 


=  0 


(12) 


(13) 


(14) 


(15) 


(16) 


second 


(17) 


(18) 


6 


1^3  +  P 


(19) 


dy 


2k 


Tl^  +  P 


(20) 


and 


(21) 


where  /S  and  yo  are  integration  constants. 

Eqns.  16  and  21  are  exact  solutions  to  the  differential 
eqns.  10  and  11.  In  fact,  boundary  conditions  are  needed  to 
evaluate  the  integration  constants  a,  /?,  Xo,  y©  as  well  as  X. 
An  important  step  in  this  procedure  is  the  recognition  of  a 
and  jS  as  functions  of  (1-M„‘)  ,  [Ref.  2],  i.e., 

a  =  ±q  (l-Ad)^  (22) 


P  =  ±C2  (1-Afi)  (23) 

The  upper  positive  sign  and  the  lower  negative  sign  in  the 
above  equations  will  be  shown  to  belong  to  Mo,&l.O  and  Me„<1.0 
respectively,  namely, 

asO  &  /JsO  for  M„al.O 

Q!<0  &  i8<0  for  M„<1.0 

See  Appendix  B. 
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The  constants,  a,  /3  and  X  can  be  shown  to  be  related  by  the 


expression 


^  - 
X 


(24) 


Since  C,  ,€2  and  X  are  positive  constants 

qCg  =  ^X  (25) 

The  usefulness  of  these  definitions  will  be  outlined  in 
more  detail  once  we  use  the  exact  solutions  from  eqns.  16  and 
21,  while  verifying  the  correct  boundary  conditions  for  the 
pressure  coefficient  Cp. 

Upon  inserting  eqns.  22  and  23  in  eqns.  16  and  21  we  have 


Factoring  out  C,  and  C2  we  obtain 


(28) 
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eqns.  32  and  33  then  follow  as 


_dl _ 


(36) 


Y  = 


/ 


dfi 


(37) 


This  set  will  be  numerically  integrated  to  obtain 
transonic  surfaces  with  their  corresponding  Cp  and  Mach  number 
profiles . 
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IV.  NUMERICAL  INTEGRATION 


We  can  rewrite  eqns.  36  and  37  in  a  more  general  form  and 
specify  the  two  integral  limits 


X 


dZ 

\/z2  ±  (1-mI)^ 


(38) 


fi 

/ 

ji 

T  (i-wi) 


dW 


(39) 


As  indicated  before,  the  upper  and  lower  sign  within  the 
integrand  root  and  the  integral  limits  is  used  for  M„sl.O  and 
M„<1.0  respectively.  In  addition,  our  flow  of  interest  is  only 
the  first  quadrant  in  X  and  Y,  which  represents  physical 
external  flows.  To  show  this  clearly,  for  the  subsonic  case 
M„<1 . 0 .  eqns.  38  and  39  reduce  to  the  form 

^  =  f  (40) 

0  -  (1-Ari)^ 


fi 

--  / 

(l-w.) 


sjw^  -  U-M^) 


_1 

fj^d-iyf)  2  (41) 
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and  for  we  have 


5 


dZ 


0  +  (1-Afi) 


2n  2 


(42) 


r  =  - 


/ 


-(!-«-) 


.2W 


^  +  (l-^») 


fj^-d-Afi)  ^  (43) 


However,  we  will  frequently  retain  the  more  general  "+" 
notation  in  our  discussions  further  on. 

Before  we  head  to  numerical  solutions,  we  should  indicate 
that  sonic  flow  (M„=1.0)  can  be  given  in  explicit  form  through 
eqns.  40,41,42,43,  i.e.,  the  resulting  equations  for  sonic 
flow  are  ^ 

(44) 


dZ 


0 


(45) 


Solving  equations  44  and  45  we  obtain  our  sonic  solution 


(46) 


(47) 
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and  inserting  eqns.  48  and  49  we  have,  [Ref.  2] 

=  1 

3  iy-Vo)^ 


(51) 


Next,  eqns.  40  to  43  are  numerically  integrated  and 
plotted  in  Figs.  1  and  2  for  various  (e.g., 
0.8, 0.9, 1.0, 1.1, 1.2)  .  Plotting  ij  vs.  Y  for  Y>0  in  Fig.  2  is  an 
appropriate  choice  for  investigating  physical  external  flows 
as  stated  earlier. 

Our  numerical  solution  has  been  based  upon  the  Newton- 
Cotes  method  of  order  4,  [Ref.  5],  which  may  be  applied  to  an 
even  number  of  subintervals,  each  described  as, 

Jfix)  dx=  7  f  (a)  +32  f|-^J+12  f|-^|+32  j+7  f  (i?) 
(52) 

(See  appendix  A) . 
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Numerical  solutions  which  express  ?  and  as  has  been 
noted,  can  be  combined  with  eqns.  30,31  and  eqns.  34,35  and 
substituted  in  eqn.  9  to  give  a  complete  solution  for  the 
perturbation  potential  in  terms  of  the  positive  constants  X  or 
(Cl,  C2)  ,  and  (Xo,  y^)  . 
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V.  BOUNDARY  SURFACES 


For  an  inviscid  flow,  the  condition  to  be  applied  at  the 
surface  of  a  solid  boundary  is  that  the  direction  of  the  flow 
velocity  vector  is  tangent  to  the  solid  surface,  [Ref.  6] .  In 
terms  of  perturbation  velocities  this  boundary  condition 


becomes 


suxface 


The  modified  velocity  perturbation  potential  as  pointed  out 


earlier  is 


thus 


<|)^  =  Wf(Y+l>-^ 


L  AC(Y+1) 


We  can  show  by  substitution  that 


m 


surface  M^(y+1) 


Now,  deriving  the  perturbed  potential  function  in  eqn.  9 
with  respect  to  y  gives 
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Eqn.  19  for  y>0  becomes 


3  ' 


(58) 


This  relation  can  be  rearranged  using  5)  and  j(3  in  eqns. 
23  to  yield 


31  and 

(59) 


Upon  using  this  equation  and  the  relation  of  eqn,  30  to  the 
derived  velocity  potential  in  eqn.  57  results  in 


3X 


2C 


1/ 


1  1  _ 


(60) 


The  last  equation  can  be  rewritten  with  the  benefit  of  the 
related  constants  Cj,  Cj  and  X  in  eqn.  25,  hence 


(61) 


Once  we  have  found  the  form  of  the  expression  for  the 
surface  boundary  condition  is  given  as 


2  1 
2  MUy+1) 


(l-A/») 


(62) 


However,  an  important  result  to  observe  is  that  ^  and  5)  in  the 
R.H.S.  of  this  equation  are  given  implicitly  as  a  function  of 
X  and  y  respectively  in  our  earlier  numerical  integration 
solution,  while  the  L.H.S.  is  expressed  by  x  and  y.  To  get 
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compatibility,  we  then  transform  the  L.H.S.  of  eqn.  62  using 


the  chain  rule 


^  ^  dX/dx 
dx  dX  dY/dy 


(63) 


Deriving  eqn.  34  with  respect  to  x  and  eqn.  35  with  respect 
y  gives 


dX  _  ^-i/3X\i 
dx  ~  ^  [  2  j 


(64) 


dY 

dy 


(65) 


therefore 


^  -  3 

dx  ~  2  dX 


(66) 


Equation  62  can  be  rewritten  using  the  above  transformation  as 


dx 


(If _ ^ 

\3/  AfJ(Y+l) 


fj3  ± 


(67) 


or 


dY 

\j  ± 


(2f _ L_ 


I  dX 


(68) 


We  are  now  in  a  position  to  compute  the  exact  boundary 
surfaces  out  of  eqn.  68.  Starting  with  the  L.H.S.  we  then 
define 


EiY)  = 

\l  ± 


(69) 
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Having  defined  eqn.  69  we  can  plot  E(Y)  vs.  Y  in  Fig.  3  for 


various  M„,  (e.g.,  0.8, 0.9, 1.0, 1.1, 1.2),  keeping  in  mind  our 

"±"  notation  for  the  supersonic  (upper  sign)  and  subsonic 
(lower  sign) . 

The  asymptotic  value  of  E(Y)  occurs  at  Yq  where 

fj  -  t(I-AL^) 


In  other  words,  the  asymptotic  value  can  be  expressed  as 
Yo=Yo(M«),  or  as  function  of  (1-M„^)  ,  which  has  been 
numerically  found  to  be 


3^0  “ 


2.4 


1-mI 


1 

2  I  6 


(71) 


Yo  vs.  Mo,  is  sketched  in  Fig.  4  to  demonstrate  an  interesting 
symmetry  with  respect  to  M„=1.0  .  Here  Yo=oo  but  note  that  for 
practical  reasons  we  can  adopt  finite  values  of  Yq  (perhaps  as 
low  as  7  to  20)  to  represent  conditions  sufficiently  close  to 
M„=1.0  . 

Defining  the  integral  of  the  R.H.S.  of  eqn.  68  as 


K{X) 


X 


(72) 


and  redoing  numerical  integration  allow  us  to  plot  K(X)  vs.  X 
for  7=1.4  in  Fig.  5.  It  is  important  to  realize  that  K(X)>0 
for  M„<1.0  and  K(X)<0  for  M„fel.O  ,  before  solving  eqn.  68. 
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Therefore,  integrating  both  sides  of  eqn.  68  yields  for  K(X)>0 


^E{Y)dY=K{X)  <73) 

and  for  K(X)<0 

y 

jE{Y)dY  =  K{X)  (74) 

i'o 

Once  we  specify  X  for  K(X)  we  can  immediately  solve  the  L.H.S. 
integrals  of  eqns.  73  and  74  for  Y,  expressing 

(^Kurface  =  (75) 


namely,  the  required  boundary  surface  for  any  chosen  M„. 

We  can  rewrite  eqns.  73  and  74  in  the  more  complete  form 
for  Mo,<1.0  as 


/ 


y 


dY 

/  fj^  -  u-m!) 


I2\^  1 


and  for 


y 


/ 


dY 

^  fj^  +  (1-M^) 


(76) 


(77) 


Eqns.  76  and  77  are  equivalent  to  eqns.  73  and  74  which 
finally  yield  the  boundary  surfaces.  These  are  shown 
nondimensionaly  in  Figs.  6  and  7  for  Mc„=0 . 8 , 0 . 9, 0 .97  and  in 
Fig.  8  for  iyi(„=l . 0, 1 . 1, 1 .2  . 


19 


The  sonic  surface  is  most  easily  verified  by  direct 
analytical  expansion  of  eqn.  77,  using  the  sonic  relations 


stated  in  eqns.  46  and  47. 


Solving  the  above  integrals  on  both  sides  gives 


(78) 


(79) 


Recall  that  Yq  for  the  sonic  case.  However,  choosing  for 
example  yo=20  gives  an  asymptotic  value  of  ?7=0.0l,  (eqn.  47), 
and  E(Yo)=1000,  (eqn.  69),  for  our  explicit  sonic  solution. 
(See  also  Figs.  2  and  3) .  Therefore,  dividing  by  a  finite  Yq 
enables  us  to  obtain  the  nondimens ional  sonic  boundary  surface 
as  a  function  of  7  . 


_r 

5^0 


1-  ^ 

(xY 

(Y^l)  \3)  1 

(80) 


This  will  be  shown  to  be  a  divergent  surface  which  increases 
the  Mach  number  to  supersonic  conditions. 
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VI .  PRESSURE  COEFFICIENT 


For  a  flow  with  small  perturbations  the 
coefficient  is  given  by,  [Ref.  8] 


A  first  order  approximation  for  the  linearized 
coefficient  in  a  two  dimensional  planar  flow  is 

^  _  2u 

^  “  U 
'*00 

Recall  the  modified  velocity  potential 

4)^.  = 


Thus 


_u  _ 

Mi(Y+l) 


and  from  eqn.  82  we  see  that 

c  = 

Mhy*l) 


Now,  deriving  the  potential  function  in  eqn 
respect  to  x  we  have 


pressure 

(81) 

pressure 

(82) 

(83) 

(84) 

(85) 

9  with 

(86) 
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Rewriting  eqn.  14  using  the  relation  for  oi  in  eqn.  22  as 


dx  \ 


(87) 


By  taking  advantage  of  |  as  stated  in  eqn.  31  and  factoring 
out  C, ,  eqn.  87  becomes 

V  V  i  (1-Mi)"  (88) 


Substituting  eqn.  88  and  the  relation  for  rj  in  eqn.  31  to  the 
velocity  potential  in  eqn.  86  yields 


4>x 


=  fi 


I  21  ) 


V  V  ±  (1-wf)^  +  (1-wf) 


(89) 


Once  again,  using  the  benefit  of  the  related  constants  Cj,  C2 
and  X  in  eqn.  25,  eqn.  89  becomes 


4)^  =  fi  7  V  ±  (1-Ari)"  +  (1-Mi) 


(90) 


therefore  eqn.  85  for  Cp  follows  as 


Cp  =  - [fi  (1-wi)'  +  (1-A/i) 


(91) 


We  must  anticipate  the  above  expression  for  Cp  to  satisfy  the 
boundary  conditions  on  external  boundary  surfaces  as  will  be 
explored  further  on. 
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VII.  THE  LOCAL  MACH  NUMBER 


The  velocity  vector  V  is  given  by  the  undisturbed  uniform 
velocity  and  the  perturbation  velocities  for  the  2D  flow  as 

V  =  £{U^+u)  +  Jv  (92) 

Thus,  we  may  obtain  the  local  Mach  number  M  in  terms  of  the 
above  velocities  and  the  local  speed  of  sound  "a". 


\V\^  _  {U„  +  u)^  + 


Neglecting  higher  powers  of  the  perturbation  velocities  gives 


= 


In  addition,  the  ratio  of  the  local  speed  of  sound  to  the 
undisturbed  uniform  speed  of  sound  becomes 


^ 

ai  "  7’o/r 


were  Tq  designates  the  stagnation  temperature. 
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Eqn.  96  can  be  replaced  using  the  well  known  relation 


l  +  UzM^ 

2 

l  +  lZlAf2 
2 


(97) 


Substituting  the  local  speed  of  sound  from  eqn.  97  to  eqn.  94 
results  as 


/ 


=  Ml 


1  + 


2u 


^1+ 

2 


(98) 


Having  stated  Cp  in  eqn.  82,  eqn.  98  is  reduced  to  the  form 


l+V^l^z  i+V^mI 
2  2 


(99) 


Rearranging  eqn.  99  and  factoring  out  yields  the  final 
expression  for  the  local  Mach  number 


1  +  M^C 


(100) 


Note  that  M=M„  when  Cp=0  as  it  should. 

This  important  relation,  combined  with  the  corresponding 
eqn.  91  for  Cp,  will  be  used  when  discussing  the  boundary 
conditions  in  the  next  chapter. 
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VIII.  BOUNDARY  CONDITIONS 


Rather  than  providing  the  boundary  conditions  for  the 
governing  differential  eqns.  10  and  11  to  evaluate  the 
constants,  we  imposed  a  solution  by  defining  a  and  /?  in  eqns. 
22  and  23  respectively.  That  is  to  say,  boundary  conditions 
have  been  already  implied  in  the  definitions  of  a  and  /5.  This 
procedure  led  us  to  reveal  the  boundary  surfaces  upon  using 
the  tangency  condition  and  compute  the  pressure  coefficient  Cp 
on  the  body  surface,  (eqn.  91) . 

Thus,  we  arrive  at  the  necessity  to  verify  the  validity  of 
the  expression  for  Cp  for  the  boundary  surface,  and  verify  the 
following  requirements: 

•  At  the  upstream  end  (X/Yo=0)  ,  we  must  require  that  Cp=0 
ensuring  the  unperturbed  Mach  number  of  the  transonic 
inflow  to  equal  M„. 

•  For  the  subsonic  inflow  at  the  downstream  end  (X  @  Y/Yo=l) 
we  must  require  that  the  local  Mach  number  saturates 
towards  M=1.0  .  This  is  because  of  the  well  known  relation 
in  gas  dynamics  for  a  convergent  flow. 

We  first  restrict  our  attention  to  the  first  requirement, 
for  M„&1.0  at  X/Yo=0 

I  =  V  =  0  (101) 

At  Y/Yo=l  eqn.  70  becomes 

fj  -  fj*  =  -(I-mJ)  2 
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Inserting  the  above  terms  of  |  and  ^  in  eqns.  90  and  91  gives 
4)^  =  -d-wi)  3  +  (1-Afi)  =  0 

and 

Cp@  (x/Y^  =  0  ,  r/y-Q  =  1)  =  0  (104) 

Therefore,  from  eqn.  100 

M®  {X/Yq  =  0  ,  Y/Y^  =  1)  =  ^„  (105) 

Hence,  eqns.  104  and  105  comply  with  the  first  requirement. 
The  complete  behavior  of  the  pressure  coefficient  and  the 
local  Mach  number  are  plotted  in  Figs.  9  and  10  for  M„=1.0  and 
M„=l.l.  In  addition,  it  is  important  to  note  that  the  boundary 
surfaces  ought  to  be  truncated  at  M„sl.2,  so  as  not  to  exceed 
the  small  perturbation  assumptions. 

Now,  for  M„<1.0  at  X/Yo=0 

^  =  0  (106) 

However,  at  Y/Yo=0 

fj  fj*  =  00  (107) 

Thus,  eqn.  91  results  as 

Cp  ®{X/Y^  =  0  ,  Y/Yq  =  0)  -  CO  (108) 

We  immediately  conclude  that  the  above  coordinate  (0,0)  can 
not  serve  as  the  upstream  start  for  the  subsonic  inflow. 
Instead,  we  must  localize  a  new  upstream  end  coordinate  along 
the  boundary  surface  whose  and  5)*  values  satisfy  eqn.  91 
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for  the  equality  Cp=0  .  Since  every  coordinate  on  the  boundary 
surface  (X/Yq/Y/Yo)  is  a  function  of  |  and  (eqn.  76)  ,  we  can 
numerically  scan  the  locus  of  all  the  coordinates,  i.e.,  the 
boundary  surface,  for  such  values  of  and  that  will 
satisfy  eqn.  91,  For  example,  scanning  the  results  for  the 
subsonic  boundary  surface  of  M„=0.8  reveals;  |=|*=- 0 . 35763  and 
^=^*=3.01393  at  the  coordinate  (X/Yo=0. 31366, Y/Yo=0. 40616)  . 
Inserting  the  above  values  of  and  rj*  in  eqn.  91  yields 
Cp=0.0003  .  Hence,  by  relocating  a  new  upstream  coordinate  for 
the  subsonic  inflow  we  are  then  able  to  satisfy  our  first 
requirement.  (This  is  shown  later  in  Figs.  11  and  12). 

We  now  consider  the  second  requirement  for  the  subsonic 
outflow  at  the  downstream  end  where 

5  -  r  = 

fj  ®(Y/Yo  =  l)  -  fj*  =  (1-m2)  2 
or  equivalently  dY/dX  ->  0 

Substituting  these  terms  for  f  and  ^  in  eqns.  90  and  91  gives 

4)^=(1-Mi)  (111) 


(109) 

(110) 


and 


-2(1-M^) 

MUy+1) 


(112) 
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Thus,  we  can  express  the  local  Mach  number  using  eqn.  100 


= 


Mt 


I  ■>  \ 

^  2{1-mI) 

1  + - 


_  (Y- 


wf  (y+1) 


(113) 


Upon  expanding  eqn. 


113  and  factoring 

.  Mi(Y-l)+2 

=  - 

m:(y-1)  +2 


out  M„^  we  obtain 

(114) 


or  equivalently 

M  =  1 


(115) 


Hence,  we  have  verified  the  second  requirement:  The  local 
Mach  number  saturates  at  the  downstream  end  to  the  value  M=1.0 
for  any  subsonic  inflow  M„<1.0  .  The  pressure  coefficient  and 
the  local  Mach  number  behavior  along  the  boundary  surface  are 
plotted  in  Figs.  11  and  12  for  Moo=0.8  and  Moo=0.9  respectively, 
showing  clearly  both  the  upstream  and  downstream  boundary 
conditions . 

Knowing  that  boundary  conditions  require  the  gradient  of 
the  potential  function  0  to  vanish  far  away  from  the  body 
gives  rise  to  an  important  question:  How  is  the  flow  pattern 
along  the  boundary  surface  propagated  into  the  flow  field?  We 
treat  this  question  by  considering  a  computationally  efficient 
flow  solver  developed  by  J.  A.  Ekaterinaris,  [Ref.  7]  .  This 
CFD  solver  consists  of  the  inviscid  Euler  equations  and  uses 
fully  implicit  two  dimensional  Crank-Nicholson  method  with  a 
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central  difference  ADI  scheme,  and  fourth  and  second  order 
numerical  dissipation,  [Ref.  9]  .  This  solver  has  been  modified 
here  to  solve  the  inviscid  flow  around  our  earlier  boundary 
surface  solution  for  the  chosen  sonic  case  M„=1.0  .  By 
creating  a  grid  mesh  which  consists  of  a  boundary  surface 
edge,  we  are  able  to  map  the  flow  pattern.  Upon  observing  the 
constant  Mach  lines  as  graphed  in  Fig.  14,  we  deduce  that  the 
flow  is  shockless  through  the  whole  transonic  range;  a  shock 
is  formed  further  downstream  at  the  local  Mach  number  M=l.9  . 
These  results  are  in  complete  accordance  with  the  inviscid 
small  perturbation  transonic  solution,  (see  Fig.  9)  .  In 
addition,  each  vector  on  Fig.  13  describes  the  resultant  flow 
direction  at  a  specified  grid  point,  (not  the  streamlines  !) , 
colored  by  Mach  levels.  Thereby,  we  may  evidence  the  far- field 
boundary  conditions 

V=C7„  u  =  v=  0  (116) 

In  other  words,  the  flow  pattern  is  propagated  out  from  the 
boundary  surface  gradually  towards  satisfying  the  far- field 
boundary  conditions  far  away  from  the  body. 

It  is  of  special  interest  to  inspect  the  dimensional 
boundary  surface  for  possible  design  applications;  For 
example,  the  dimensional,  subsonic,  boundary  surface  of 
M5„=0.8  is  obtained  by  multiplying  the  nondimensional 
coordinates  (X/Yo(o.8)/ Y/Yoto.s))  .84  ,  i . e. ,  the 
dimensional,  downstream  end  coordinate  results  as  (Xi,2.84), 
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as  plotted  in  Fig.  15a.  Similarly,  multiplying  the 
nondimensional,  sonic  coordinates  (X/Yo(i.o),y/Yo(,.o))  by  Yo(,.o)=14, 
yields  the  dimensional,  sonic,  boundary  surface,  where  the 
dimensional,  upstream,  starting  coordinate  becomes  (0,14). 
Now,  translating  the  dimensional,  sonic,  boundary  surface, 
(Pig.  15b),  by  (AX,AY)  =  (X,,  -  (14-2 . 84) )  ,  to  match  the 
downstream  end  coordinate  of  the  above  dimensional,  subsonic, 
boundary  surface,  "patches"  both  dimensional  boundary 
surfaces,  to  form  a  2D  planar  upper  surface  as  shown  in  Fig. 
15.  We  may  use  this  technique  to  "patch"  any  subsonic  boundary 
surface  in  order  to  form  our  transonic  upper  surface  of 
interest . 
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IX.  CONCLUSIONS 


We  have  shown  that  an  exact  solution  for  the  transonic 
equation  satisfies  the  complete  behavior  of  the  pressure 
coefficient,  and  the  local  Mach  along  a  2D- surface  for  planar, 
external  flow.  It  is  of  special  interest  to  inspect  the 
dimensional  boundary  surface  for  possible  design  applications; 
For  example,  the  dimensional,  subsonic,  and  sonic  boundary 
surfaces  are  obtained  by  multiplying  the  nondimens ional 
coordinates  (X/Yo,Y/Yo)  ,  by  the  corresponding  Yo.  Translation 
of  the  dimensional,  sonic,  boundary  surface,  to  match  the 
downstream  end  coordinate  of  the  above  dimensional,  subsonic, 
boundary  surface,  "patches"  both  dimensional  boundary 
surfaces,  to  form  a  2D  planar  upper  surface  as  shown  in  Fig. 
15.  This  technique  may  serve  as  a  powerful  tool  for  the 
implementation  of  supercritical,  shockless,  planar  surface 
designs,  because  of  the  saturation  of  yj  with  Y  (Fig.  2),  we 
can  also  describe  certain  internal  flows  such  as  planar 
converging,  diverging  tunnel  walls.  We  may  also  extend  these 
results  to  axisymetric  transonic  flow. 
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APPENDIX  A 
COMPUTER  PROGRAMS 

Q* ********************************************************** 

*  THIS  FORTRAN  PROGRAM  USES  NUMERICAL  INTEGRATION  BASED  ON 
THE  NEWTON- COTES  FORMULA  TO  SOLVE  NUMERICALLY  EQNS.  40  AND  42. 

THE  OUTPUT  IS  IN  THE  FORM  OF  f  vs.  X 


program  KSI 
PARAMETER  (N=5) 

IMPLICIT  REAL*8 (A-H,0-Z) 

IMPLICIT  REAL* 8 'M) 

IMPLICIT  REAL* 8 (K) 

INTEGER  INDEX, I 

DIMENSION  X (N) , XDOT (N) , SAVEX (N) , SAVED (N) , KSI (N) 

DIMENSION  MACH (N) , XI (N) , X2 (N) 

MACH(l) =0.8D0 
MACH(2)=0.9D0 
MACH(3)=1.0D0 
MACH (4 ) =1. IDO 
MACH(5) =1,2D0 

KSI(1)=0,0D0 
KSI (2)=0.0D0 
KSI (3) =1.0D-3 
KSI(4)=0.0D0 
KSI(5)=0.0D0 

Q* ************************************************* 

C  MAIN  PROGRAM 

Q*  ************************************************  * 

DO  100  1=1,5 
IF(I.LT.3)  THEN 
H=-0.001D0 
GOTO  5 
END  IF 
H=0.01D0 

5  IF(I.LT.3)  THEN 

TF=- (l.ODO-MACH(I) **2) 

GOTO  10 
END  IF 
TF=120.0D0 

10  IF(I.LT.3.AND,KSI(I) .GT.TF.OR. 

I.GE.3 .AND.KSI (I) .LE.TF)  THEN 
30  do  50  index=l,5 

CALL  DERIVative (KSI, N, I, XI, X2, XDOT, MACH) 

CALL  NEWTONCOTES (KSI, XDOT , X , SAVEX , SAVED , H , N , INDEX , I ) 
50  end  do 
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WRITE(10+i,*)  X{I),KSI(I) 

GOTO  10 
END  IF 
100  END  DO 
STOP 
END 

c* *************************** ******************************* 

c  SUBROUTINE  DERIVative  * 

Q*********************************************************** 

SUBROUTINE  DERIVative (KSI , N, I , XI , X2 , XDOT, MACH) 

IMPLICIT  REAL*8 (A-H,0-Z) 

IMPLICIT  REAL* 8 (M) 

IMPLICIT  REAL*8(K) 

DIMENSION  XDOT (N) , XI (N) , X2 (N) , KSI (N) , MACH (N) 

GO  TO  (20,20,10,10,10)  ,I 
10  XI (I) =KSI (I) **2+ (l.ODO-MACH(I) **2) **2 
X2 (I)=Xl(I)**{1.0D0/3.0DO) 

XDOT(I)=1.0D0/X2 (I) 

GOTO  30 

20  X1(I)=- (KSI(I)**2- (1.0D0-MACH(I)**2)**2) 

X2 (I)=-X1(I) ** (1.0D0/3.0D0) 

XDOT(I)=1.0D0/X2 (I) 

30  RETURN 
END 

(3*********************************************************** 

c  SUBROUTINE  NEWTONCOTES  * 

C*** ****** ************************************************** 

SUBROUTINE  NEVJTONCOTES  ( KSI ,  XDOT ,  X ,  SAVEX ,  SAVED ,  H ,  N ,  INDEX ,  I ) 
IMPLICIT  REAL*8 (A-H,0-Z) 

IMPLICIT  REAL* 8 (M) 

IMPLICIT  REAL* 8 (K) 

DIMENSION  X (N) , XDOT (N) , SAVEX (N) , SAVED (N) , KSI (N) 

GO  TO  ( 1 , 2 , 3 , 4 , 5 )  , index 

1  SAVEX (I) =X (I) 

SAVED(I)=XDOT(I) 

KSI (I) =KSI (I) +0.25D0*H 
RETURN 

2  SAVED ( I ) =7 . ODO* SAVED ( I ) +32 . ODO*XDOT ( I ) 

KSI (I) =KSI (I) +0.25D0*H 

RETURN 

3  SAVED(I)=SAVED(I) +12 . ODO*XDOT(I) 

KSI (I) =KSI (I) +0.25D0*H 

RETURN 

4  SAVED ( I ) =SAVED ( I ) +32 . ODO*XDOT ( I ) 

KSI (I) =KSI (I) +0.25D0*H 

RETURN 

5  X(I)=SAVEX(I) + (H/90.0D0) * (SAVED (I) +7 . ODO*XDOT ( I) ) 

RETURN 

END 

Q*********************************************************** 
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************************************************************ 
THIS  FORTRAN  PROGRAM  USES  NUMERICAL  INTEGRATION  BASED  ON 
THE  NEWTON-COTES  FORMULA  TO  SOLVE  NUMERICALLY  EQNS .  41  AND  43. 

THE  OUTPUT  IS  IN  THE  FORM  OF  ^  vs .  Y 
************************************************************ 

program  ETTA 
PARAMETER  (N=5) 

PARAMETER  (K=2000) 

IMPLICIT  REAL*8 (A-H,0-Z) 

IMPLICIT  REAL* 8 (M) 

INTEGER  INDEX, I, J 

DIMENSION  Y(N) ,YDOT(N) ,SAVEY(N) , SAVED (N) ,TF(N) ,MACH(N) 

DIMENSION  ETA(N)  ,ETA11(K)  ,Y11{K)  ,Y1(N)  ,Y2{N) 

MACH(1)=0.8D0 

MACH(2)=0.9D0 

MACH(3)=1.0D0 

MACH(4)=1.1D0 

MACH(5)=1.2D0 

H=0.01D0 

c********* ************************************************** 

C  INTEGRAL  FINAL  LIMITS  * 

c* ********************************************************** 

DO  7  1=1,5 

IF(MACH(I) .EQ.IODO)  THEN 
TF(I)=0.0D0 
GOTO  7 
END  IF 

IF(MACH(I) .GT.l.ODO)  THEN 

TF(I)=(- (1.0D0-MACH(I)**2) ) ** (1 . ODO/3 . ODD) 

GOTO  7 
END  IF 

TF(I)=(1,0D0-MACH(I)**2)** (1 . ODO/3 . ODO) 

7  END  DO 

c* ********************************************************** 
C  MAIN  PROGRAM  * 

c********* **************  ************************************ 

DO  100  1=1,5 
J=1 

ETA (I) =14. ODO 

Y (I) =2. ODO/ (ETA (I)** (1 . ODO/2 . ODO) ) 

10  IF(ETA(I) -H.GT.TF(I) )  THEN 
20  do  50  index=l,5 

CALL  DERIVative ( ETA, N, I , Y1 , Y2 , YDOT, MACH) 

CALL  NEWTONCOTES (ETA, YDOT, Y, SAVEY, SAVED, H, N, INDEX, I) 
50  end  do 

ETA11(J)=ETA(I) 

Y11(J)=Y(I) 

J=J+1 
GOTO  10 
END  IF 
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H=0.001DO 

IF(ETA(I) -H.GT.TF(I) )  THEN 
GOTO  20 
END  IF 

80  DO  90  L=J-1, 1, -1 

WRITE (20+i,*)  Yll (L) ,ETA11 (L) 

90  END  DO 
100  END  DO 
STOP 
END 

C* ********************************************************** 

C  SUBROUTINE  DERIVative  * 

Q* *******  **************************************  ************* 

SUBROUTINE  DERIVative (ETA, N, I , Y1 , Y2 , YDOT, MACH) 

IMPLICIT  REAL*8 (A-H,0-Z) 

IMPLICIT  REAL*8(M) 

DIMENSION  YDOT(N) ,Y1(N) ,Y2 (N) ,ETA(N) ,MACH(N) 

GO  TO  (20,20,10,10,10)  ,I 
10  Y1 (I)=ETA(I) **3+ (l.ODO-MACH(I) **2) 

Y2 (I)=Y1(I)** (1.0D0/2.0D0) 

YDOT(I)=1.0D0/Y2 (I) 

GOTO  30 

20  Y1(I)=ETA(I)**3- (l.ODO-MACH(I) **2) 

Y2 (I)=Y1(I)** (1.0D0/2.0D0) 

YDOT(I)=1.0D0/Y2 (I) 

30  RETURN 
END 

Q* ********************************************************** 

C  SUBROUTINE  NEWTONCOTES  * 

C* ********************************************************** 

SUBROUTINE  NEWTONCOTES ( ETA , YDOT , Y , S AVE Y , SAVED , H , N , INDEX , I ) 
IMPLICIT  REAL*8 (A-H,0-Z) 

IMPLICIT  REAL* 8 (M) 

DIMENSION  Y (N) , YDOT (N) , SAVEY (N) , SAVED (N) , ETA (N) 

GO  TO  ( 1 , 2 , 3 , 4 , 5 )  , index 

1  SAVEY (I) =Y (I) 

SAVED(I)=YDOT(I) 

ETA(I) =ETA(I) -0.25D0*H 
RETURN 

2  SAVED(I)=7.0D0*SAVED(I) +32 . ODO*YDOT (I) 

ETA(I) =ETA(I) -0.25D0*H 

RETURN 

3  SAVED ( I ) =SAVED ( I ) +12 . ODO* YDOT ( I ) 

ETA(I) =ETA(I) -0.25D0*H 

RETURN 

4  SAVED(I)=SAVED(I) +32 . ODO*YDOT(I) 

ETA(I)=ETA(I) -0.25D0*H 

RETURN 

5  Y(I)=SAVEY(I) +(H/90.0D0) * (SAVED (I) +7 . ODO*YDOT (I) ) 

RETURN 

END 
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APPENDIX  B 
FORM  OF  a  AND  /3 


The  expressions  for  a  and  j8  as  introduced  in  [REF.  2]  are 

a=±q(l-^2)"  (117) 


p  =  +C2  (1-wf)  (118) 

for  0.8sM„sl.2  . 

Upon  inserting  eqns.  117  and  118  in  eqns.  16  and  21  yields 

X=  [-—====  (119) 

0 


(120) 


These  relations  are  numerically  integrated  and  plotted  in 
Figs.  16  and  17.  We  observe  the  same  solutions  for  M„sl.O  as 
expected,  but  distinguish  a  different  set  of  solutions  for  the 
subsonic  case  M„<1.0  . 

Applying  the  same  procedure  as  in  Chap.  IV,  we  easily 
compute  the  expression  for  the  boundary  surfaces 


(121) 
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The  numerical  integration  of  eqn.  121  is  shown  in  Fig.  18  for 
various  (e.g.,  0.8,0.9,0.97,1.0,1.1,1.2).  The  diverging 
contours  of  the  subsonic  solutions  appear  to  give 
contradictory  results  once  investigating  the  pressure 
coefficient  behavior 


C-p  = 


-2 


+  {I-mD"  +  (l-AC) 


(122) 


along  the  boundary  surface.  This  implies  as  shown  in  Fig.  19 
that  the  boundary  surface  within  the  portion  of  Cp>1.0  is  not 
a  realistic  solution. 
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APPENDIX  C 
FIGURES 


Figure  1.  Numerical  integration  of  eqns.  40  and  42. 
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ETTA(Y)forM=0.8-1.2 


Figure  2.  Numerical  integration  of  eqns.  41  and  43. 
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E(Y)  vs.  Y  for  M  =0.8- 1.2 


Figure  3.  Numerical  solution  of  E(y) . 


Figure  4.  Equation  71. 


K(X)vs.X  for  M„=0.8-1.2 
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Figure  5.  Numerical  integration  of  eqn.  72. 
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Figure  6.  Numerical  integration  of  eqn.  76. 
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Figure  7.  Subsonic  boundary  surfaces. 


Figure  8.  Numerical  solution  of  eqn.  77 
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Figure  9.  Cp  and  local  Mach  for  Ma,=1.0  . 


Figure  10.  Cp  and  local  Mach  for  Moo=l.l  . 
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Figure  11.  C  and  local  Mach  for 


Figure  12.  C  and  local  Mach  for 


v  ,2  S.*'  «■  p  oi- 

,  2  <fo!::,'P  o 


i  i  i  i  i  i  i 

-S  <  >  >  i  J  f 


f  I  i  '  ^ 

/  I  i  f  r 


/f??  ip-  /  i 


f  ?  ;  /  ^  fQ^  J  J  f  I  1  d 

I  f  i  •?  I  H  f  ! 


/  I  I  /  f  //////  /  i 

I  I  I  I  J '  I  I  I  I  I  »|  I”  if  ^|■'  ,//  >  / 

i  i  /  /  I  M  I  ^  '  f,  / 

i  /  (W  /  /  -v  f  ^  f  ^/f  :  /  ' 

i  i  i  ii  I  ^  ^  i  I  i  f  h  .:  j/ 


f  'f  '1 


4/  i'//  . 

fj 


«;  |;  Ji  ^  iV  4^  ^  «  1=  ^  ;f  '  /'  4  i;  X:  c*.’  4?  ^ 

iJhiUiiffffHHiifff 

'll  ?  I  f  i’  ?  f  f  f  i| ,!' / 

1 1  j  n H n^ahiiiii 

§  i  i  i  h  f  u  n  f I  f  '(hiii^ 

I  .'  j  :  J  is!#  t*  i  i'  s  j  I  *  ’  f  / 

5  ^  i  M  f  f  i'  r  =■  f  ;  i  j  .!  «  f  ^  i 

6  I  j  I  I  H  M  H  M  ^  ;'  /  I  I  /  /  / 
HI*..'  j!  -i  il  i'  i  <•■  I  i  I  r 

if  i  l  i  i  n  i  fH(fi'l!ljlj 

I n  '  i  >.  n  I  1 1  If  f  I f f f f I r 

^  i  I  J  i  J  ^  f  li  *  >  ^  H  n  !  :  H 

?  i  i  i  i  I  i  i  s  i  ■<  M  1 1 ;  1 1  f  M 

II  i  I  1  I  I  H  ^  U  H  ^  M  f  I!  t 
I  I  I  i  i  j  i  i  1  i  I  I  I  M  M  H  f  I 
|.  I  I,  I  I  I  I  i  i  I  I  I  I  f  •>  i  1 1  -M 

;  !  :  I  I  :  i  i  1  •  !  j  1  ■  i 

I  I  4  j  ?• 

i  I  5  ■  !  i  i  f  I  If-  if  i  I  ' 

fj  I  I  ?,  I,  J  ^  I  I  I  I  I  I  I  i  I  i  I  I 

if  if  i  '  I  i  1  f  :  f  ■  i 

ff  ii  X  i  4  J  I  4  ^  i  I  I  I  I  I  I  I  1  I  I  I 

J  ;  ;  I  s  *  ?  I  j  i  :  f  j  I  5  : 

I  f|  ?  fj  fj  ’‘i  ‘i  U  'k  ^  t  I  a  I  I  I  I  I  I  i 

u  I  I  4  jj  i  i|  ^  f j  ■!  I  ^  I  I  I  j  I  I  I 


|.  .j  fj  i|.  j 


1  a 

f#  h  i  u 

t 

4  f.' 

1  i  1 

: 

U  1 

.««L 

i  i  i  i  i 

■-!•  . 

A  ^  *1 

luj'i 

H  ^  ^ 

^  H  M  1 

5 

i>  i 

•  i 

i  1  ?  ' 

•  : 

iUji 

j  i  St 

C£ 

I  §88; 

2888; 

O  *  ' 

Q  o  O  C  ■: 


8  0  ^  ^  X  ^  > 

S  ?5  ^  X  j«;  v<  >k:  > 

'■'''  o  “"  <  *  <■ 


?,S88S88S 


'''*f  •■  ^  S/»  Its  4"  >. 

?.?  ■^^  O  •'*<  ‘•**  •••*  '**  •••■ 


«*■  “5  »«  **•'  ‘C-’  ^2*  CJ?  O'  ‘h¥’' 

-*•  .*,*  O  fi'i  xn  nCT';  W  9^ 

....  6i  ?>]  H  c4  fs  ^  ri 


Figure  13 ..  Velocity  yectores  for 
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Figure  14.  Mach  lines  for  M„=1.0 


Figure  15.  Upper  transonic  surface  for 


Figure  15a.  Upper  subsonic  surface  for 
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Figure  15b.  Upper  supersonic  surface  for 
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Figure  16.  Numerical  integration  of  eqn.  119. 
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Figure  17.  Numerical  integration  of  eqn.  120. 
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Figure  18.  Numerical  integration  of  egn.  121. 
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